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Abstract We study numerically nonuniform quantum turbulence of coflow in 
a square channel by the vortex filament model. Coflow means that superfluid 
velocity Vg and normal fluid velocity flow in the same direction. Quantum 
turbulence for thermal counterflow has been long studied theoretically and 
experimentally. In recent years, experiments of coflow are performed to observe 
different features from thermal counterflow. By supposing that Vg is uniform 
and Vn takes the Hagen-Poiseiulle profile, our simulation finds that quantized 
vortices are distributed inhomogeneously. Vortices like to accumulate on the 
surface of a cylinder with Vg ~ d„. Consequently, the vortex configuration 
becomes degenerate from three-dimensional to two-dimensional. 
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1 Introduction 

Quantum turbulence is one of the most important issues in low temperature 
physics and has been studied theoretically and experimentally for thermal 
counterflow where superfluid and normal fluid flow oppositely. In recent years, 
experiments of coflow where superfluid and normal fluid flow in the same 
direction are performed to observe different features from thermal counterflow 
[T]. For example, the vortex line density is proportional to the 3/2 power of 
the velocity and independent of temperature. The motivation of the present 
paper is to find numerically some behavior of vortices characteristic of coflow. 
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According to the two-fluid model, superfluid ^He at a finite temperature 
consists of an intimate mixture of a viscous normal fluid component (density 
Pn and velocity t;„) and an inviscid superfluid components (density ps and 
velocity v^) [!]■ In superfluid ^He the circulation of a superfluid vortex, called 
a quantized vortex, is quantized by the quantum circulation n = h/rrii, where 
h is Planck’s constant and m 4 is the mass of a ^He atom. Quantum turbulence 
generally occurs by tangling of quantized vortices [3] . 

We perform numerical simulation for coflow in a square channel; super¬ 
fluid velocity is prescribed to be uniform flow and normal fluid velocity to 
be Hagen-Poiseuille flow. Vortices are distributed inhomogeneously and like 
to accumulate through the mutual friction on the surface of a cylinder where 
superfluid velocity equals normal fluid velocity. As a result, the vortex config¬ 
uration becomes degenerate from three-dimensional to two-dimensional. How 
strongly the vortices accumulate depends on temperature and the averaged 
velocity. 

The contents of this paper are as follows. In Sec. [2l we shall clarify the 
formulation of the model and introduce the equation of motion. In Sec. O we 
show the results for coflow. Section 0] is devoted to the conclusion and the 
future work. 


2 Formulation 


We perform numerical simulation by using the vortex filament model with the 
full Biot-Savart law [4115] . A point s on the vortex filament is represented in 
a parametric form s = where t is time and ^ is the one-dimensional 

coordinate along the filament. The equation of motion of s is given by 

s = Us -I- as' X (u„ - Vs) - a's' x [s' x (u„ - ««)], (1) 


where a and a' are the temperature-dependent mutual friction coefficients, 
and the prime denotes derivatives of s with respect to 

The first term of the right-hand side of Eq.(l) is superfluid velocity, given 

by 

'^s — '^s,uj '^s,b H" (2) 

Here Vg^ui is the velocity field caused by vortex filaments, Vs^b the boundary 
induced field, Vg^a the applied uniform velocity field. The velocity field is 
represented by the Biot-Savart law 


'^s,uj{r) 


K f {si — r) X dsi 


( 3 ) 


where Si refers to a point on the vortex filament and the integration is per¬ 
formed along the vortex filaments. The velocity field Vg^^ is obtained by a 
simple procedure; it is just the field produced by an image vortex that is con¬ 
structed by reflecting the filament into the surface and reversing its direction. 
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The second and third terms of the right-hand side of Eq.(l) are caused 
by mutual friction. The second term makes a vortex balloon out or collapse 
inward. As discussed in [6], when the relative velocity Vns = Vn — hows 
against the mutual friction always shrinks the vortex line locally. On 

the other hand, Vns flowing along Vg^ui yields a critical radius of curvature Re- 
When the local radius i? at a point on the vortex is smaller than Rc, the vortex 
shrinks locally, while the vortex balloons out when R > Rc- 

We prescribe the Hagen-Poiseuille profile Up for though the normal 
fluid flow in the experiment [I] may be actually turbulent. When the normal 
fluid flows along the x direction, the x component of Vn is given by 

OO 

Up{y,z) = uo 

m—1,3,5,... 

where uq is a normalization factor and a and b are halves of the channel width 
along the y and z axes, respectively [7]. 

To characterize the development of vortices, we introduce the vortex line 
density (VLB) as 

where the integral is performed along all vortices in the sample volume 17. 

In this study, our calculation is performed under the following conditions. 
A vortex filament is represented by a string of discreet points. The numerical 
space resolution, namely the minimum distance between neighboring points, is 
= 8.0 X 10“"^ cm, the time resolution is Z\ = 1.0 x 10“^ s, and the computing 
box is 0.1 X 0.1 X 0.1 cm^. We consider the case of Vn = Vg, where ir„ and Vg 
are the spatially averaged normal fluid velocity and the spacial average of the 
applied superfluid velocity I’s.a, respectively. The periodic boundary conditions 
are used along the flow direction x, whereas the solid boundary conditions are 
applied to channel walls. The effects of reconnection is artificially performed, 
whenever two vortices approach more closely than A^. The initial state consists 
of eight randomly oriented vortex rings of radius 0.023 cm (Figl3[a]). 


1 - 


cosh(m7rz/2a) 


cosh(m7r6/2a) 


cos{rmTy/2a) 


(4) 


3 Results and Discussion 

3.1 Distribution of localized vortices and development of vortex 

We perform numerical simulation for coflow in a square channel and find two 
different states of vortices, namely the diffusive state and the localized state. 
The features of two states are shown in the figures; FiglT] shows the time 
development of the VLD and Fig 15] shows the snapshots of the vortices. In 
the diffusive state, the VLD oscillates irregularly and the vortices are diffusive 
(FigUfa] and Figl5[a]). In the localized state, the VLD just increases (FigUb]) 
and the vortices localize in a region shown in Figl51b],[c]. 
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Fig. 1 Temporal development of the VLD shows two different kinds of behavior depending 
on the vortex states. In [a] with T=1.35K and v'n(= Vs) = 1.0 cm/s the VLD oscillates 
irregularly, when the vortices are diffusive as shown in Fig[2[a]. In [b] with T=1.95 K and 
Vn = 1.0 cm/s the VLD just increases, when the vortices are localized as shown in Figs[2Ib] 
and [c]. (color online) 



Fig. 2 Snapshots of vortices viewed along the flow direction; [a] is a snapshot at t = 15 
s of the dynamics of Fig^a], [b] is a snapshot at t = 20 s of the dynamics of Fig^b]. 
The distribution of vortex has two states, which are diffusive [a] and localized [b]. [c] is the 
snapshot of [b] viewed from the side. 


In the localized state, the development of the vortices consists of two stages. 
As shown in FiglTJb], the first stage is 0 s < t < 0.5 s and the second is 0.5 s < t. 
In the first stage (Fig[3Ja],[b],[c]), the vortices repeats lots of reconnections, 
while they are attracted to a localized region by the mutual friction. Some 
small vortex loops are made by reconnections and balloon out to a localized 
region. Then, the VLD increases rapidly. In the second (Figl3[d]), the vortices 
protruding from the localized region towards the walls continue to extend and 
wrap the localized region. Since the reconnections seldom occur in this stage, 
the VLD increases much slower than that of the first stage. 

The vortex distribution of coflow is different from that of nonuniform ther¬ 
mal counterflow [8]. This is because that the coflow has the region where the 
mutual friction dose not work. To see where the vortices accumulate, we con¬ 
sider the localized induction approximation (LIA). Then, the second term of 
the right-hand side of Eq.(I) is as' x (u„ — Vg^a — Ps' x s"), where P is the 
quantity proportional to the quantum circulation . The mutual friction does 
not work in the region where Vn — Ws,a — Ps' x s" vanishes. For the sake of 
















Quantum turbulence in cofiow of superfluid ^He 


5 






Fig. 3 Simulations of 
the time development 
of vortex tangle in 
coflow viewed along 
the flow direction (T 
= 1.95 K, Vn = 1.0 
cm/s); [a] £ = 0 s, [b] 
£ = 0.1 s, [c] £ = 0.2 s, 
[d] £ = 0.75 s 


simplicity, if the term jSs' x s" is negligible, the region is like a square pipe 
as shown by the square dots in FigU) However, this region is modified to the 
round dots by the term fis' x s"; the position around a corner is shifted in¬ 
ward because the radius of curvature is small and the position around a side 
is shifted outward because that is large. 
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Fig. 4 The region where mutual friction does not work (color online) 
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3.2 Parameter dependence of localization 


Whether the vortices are localized or diffusive depends on the parameters. We 
introduce a dimensionless variable ^in/iout to characterize the behavior of the 
vortices; 


Lin ^ /Cin 
Lout Jr . 

^out 


( 6 ) 


Here Lin is obtained by the integration along all vortices Lin in the cylindrical 
region of the radius 0.045 cm whose central axis is along that of the pipe and, 
Lout is obtained by the integration along all other vortices Lout- The time 
development of Lin/Lout is shown in Figl5[a]. We define that the vortices are 
localized when Lin/Lout exceeds 5 after the time passes sufficiently, and that 
the vortices are diffusive when Lin/Lout is less than 5. The value of Ljn/Lout 
depends strongly on temperature and the averaged velocity v{= Vn = Vg), 
and the resulting phase diagram is shown in Figl^b]. The vortices like to be 
localized for higher temperature and faster velocity, and diffusive for lower 
temperature and slower velocity. 

Mutual friction depends on temperature and velocity, because the coeffi¬ 
cient a is dependent on temperature and increasing the velocity makes the rel¬ 
ative velocity faster everywhere. If temperature is lower and velocity is slower, 
the velocity caused by mutual friction is dominated by the self-induced veloc¬ 
ity /3s' X s" and the vortices almost move freely. Consequently, whether the 
vortices are localized or diffusive depends on temperature and the velocity. 
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Fig. 5 [a]: the time development of dimensionless quantity Lin/Lout, [t>]- its dependence 
on temperature and the velocity. In [a], the round dots are for T = 1.95 K and Vn = 0.6 
cm/s, the square dots are for T = 1.35 K and iTn = 1.0 cm/s (color online) 


4 Conclusions 

We found that the vortex development of coflow has two states; one is diffusive 
and the other is localized. At the localized state, the vortices accumulate on 
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the surface of cylinder where the mutual friction does not work. Whether the 
vortices are localized or diffusive depends on temperature and the averaged 
velocity. 

The future work is to perform a numerical simulation of coflow under tur¬ 
bulent normal fluid flow. In the experiment [T], the normal fluid flow seems 
turbulent. We will reproduce this situation and confirm this result by calcu¬ 
lating the motion of the vortex filament with the turbulent normal fluid flow. 
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